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Abstract 

This is the first of two lectures on the technique of dimensional continuation employed by Brown, 
Preston, and Singleton (BPS) to calculate such quantities as the charged particle stopping power 
and the temperature equilibration rate in a plasma. In this exposition we will examine some of 
the more basic points of dimensional continuation, with an emphasis on the Coulomb logarithm for 
electron-ion temperature equilibration. Dimensional continuation, or dimensional regularization as 
it is more properly known in quantum field theory, was originally developed as part of the renormal- 
ization procedure for the theories of the electroweak and other fundamental interactions in particle 
physics. Dimensional continuation is so general, in fact, that any theory can be unambiguously 
lifted to dimensions beyond three, and therefore the technique is powerful enough to apply in many 
other settings. The technique, however, is not well known outside the field theory and particle 
physics communities. This exposition will therefore be self-contained, intended for those who are 
not specialists in quantum field theory, and I will either derive or motivate any requisite field the- 
ory results or concepts. Of particular relevance is the analogy between the Coulomb logarithm as 
calculated by Lyman Spitzer on the one hand, and the Lamb shift as calculated by Hans Bethe on 
the other. While dimensional continuation is a well developed and a thoroughly tested method for 
regularizing any quantum field theory, BPS employs the method in a novel way that provides the 
leading and subleading behavior for processes that involve competing disparate energy or length 
scales. In particular, BPS calculated the temperature equilibration rate to leading and next-to- 
leading order in the plasma number density for any two species in a plasma that are in thermal 
equilibrium with themselves, but not necessarily with each other. No restriction is made on the 
charge, mass, or the temperature of the plasma species. It is, however, assumed that the plasma is 
not strongly coupled in the sense that the dimensionless plasma coupling parameter g = e 2 K D /AitT 
is small, where k d is the Debye wave number of the plasma. To leading and next-to-leading order 
in this coupling, the temperature equilibration rate is of the generic form T = Ag 2 ln{Cg}. The 
precise numerical coefficient A in front of the logarithm has been known for some time, while BPS 
have recently computed the constant C under the logarithm. It should be emphasized that the 
BPS result is not a model, but rather it is an exact calculation of the leading terms in a well- 
defined perturbation theory. This exact result differs from approximations and models given in the 
literature. 
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I. INTRODUCTION AND CONTEXT 



This is the first of two lectures on a new technique for calculating the temperature equi- 
libration rate between electrons and ions in a weakly to moderately coupled fully-ionized 
plasma, exact to leading and next-to-leading order in the plasma number density. This cal- 
culation was first performed in Section 12 of Ref. [lj], a work whose primary focus was the 
charged particle stopping power in a plasma. 1 This paper assumed familiarity with a num- 
ber of field theory concepts, and Section 12 relied heavily on the charged particle stopping 
power results derived in previous sections of that work. In contrast, these lectures will be 
self-contained. I will either derive or motivate the requisite field theory background for a 
complete reading of Ref. [l], with an emphasis on some of the subtleties of the calculational 
techniques and the concepts behind them (this lecture). More to the point, since the rate 
calculation stands on its own, it should be presented on its own (the following lecture). In- 
deed, since the calculation of the rate is somewhat less involved than that of the stopping 
power, it more clearly illustrates the tools and concepts imported from field theory. 

In addition to clarifying the method of dimensional continuation, this lecture will place 
Ref. [l| in the context of more familiar and traditional approaches to the rate problem. 
In particular, I will show that dimensional continuation can be viewed as a systematic 
implementation of the approach based on convergent kinetic equations. Finally, in in the next 
lecture, I will go on to derive the main result from Section 12 of Ref. [l], the rate coefficient 
(II. 3D of this lecture. By working in the Born approximation, and adopting the methods 
of Ref. jsj], I will derive this result in a much simpler manner than originally presented in 
Ref. [lj]. While Lectures I and II are self-contained, they are complementary and should be 
read as a unit. 2 

The strategy employed by Ref. [l|, hereafter referred to as BPS, consists of two steps. 
First, we will find a dimensionless parameter in which to perform a controlled perturbative 
expansion of the rate, expanding to leading and next-to-leading order in this parameter. 



1 For a short cursory version of this work, see Ref. . 

2 In other words, Lecture I (this lecture) consists of the basic theory and techniques behind dimensional con- 
tinuation, while Lecture II (next lecture) will be a specific calculation in the extreme quantum limit: only 
by performing a calculation can one understand the underlying ideas of the calculation. Ideally, I would 
then like to continue these first two lectures with three additional ones. In Lecture III, I would present the 
full calculation of the temperature relaxation rate performed in Ref. [lj], valid beyond the Born approxi- 
mation. The full calculation is accurate for a weakly to moderately coupled plasma in both the classical 
and quantum limits, and any regime in-between, regardless of the mass and temperature difference of the 
plasma species. I would then present the details of the full quantum corrections in Lecture IV. Finally, in 
Lecture V, I would come full circle, simplifying the general equilibration rate calculated in Lectures III and 
IV to obtain the extreme quantum limit of Lecture II, equation (|1.3|) of this lecture. This would provide 
two independent calculations of the rate (|1.3|) . but alas, time does not permit these last three supplements. 
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Second, we will deploy a technique from quantum field theory that will allow us to calculate 
the coefficients of these leading and subleading terms exactly. The exact leading order term 
is not very difficult to find, and has been known since the classic work of Spitzer. The 
next-to-leading order term, on the other hand, was not known exactly until the recent BPS 
calculation. The third-order term provides an estimate of the error of the calculation. When 
the plasma is weakly to moderately coupled, the error will be small and the rate will be 
approximated quite accurately by the first two terms of this expansion. 

To calculate the expansion coefficients, we will exploit a field theory technique known 
as dimensional regularization (or dimensional continuation, as I will call it here). This 
application of dimensional continuation is quite different from its intended purpose in the 
renormalization procedure of quantum field theory. Dimensional continuation was originally 
developed as an elegant regularization scheme in which the fundamental symmetries of a 
field theory could be maintained while still rendering finite the otherwise infinite integrals 
that arise when calculating Feynman diagrams. I will show how this technique can be used 
in a novel fashion to extract next-to-leading order physics that has, until now, remained 
inaccessible. In other words, I will show how dimensional continuation provides an exact 
result for the corresponding Coulomb logarithm of the process in question. I will also take 
the opportunity to correct a small algebra mistake for the electron-ion equilibration rate 
presented in Section 12 of Ref. 

A. The Problem 

The general formalism starts with a plasma composed of multiple species labeled by 
an index b, the various species being delineated by of a common electric charge e& and a 
common mass m&. 3 Each species is assumed to be in thermal equilibrium with itself at 
temperature with a spatially uniform number density n&. I will drop the subscript on 
the charge of the electron and write e e i ect ron = ~ e (with e > 0), although the electron mass 
will be denoted by m e , and the number density and temperature of the electron plasma 
component by n e and T e respectively. I will use a lower case subscript i to denote a single 
ion species of charge = Z{ e, mass raj, number density rij, and temperature Tj. I will 

3 The final constraint of BPS is that the ionization state of each component species does not change. In other 
words, the charges can be expressed as e& = Z^ e with Z\> fixed. While this simplifying assumption has its 
limitations, it facilitates the analytic calculation of the stopping power and temperature equilibration rate. 
However, one can take the Z\, to be fractional to mock-up ionization in a simple manner. Or better yet, 
it should be possible to combine the BPS results with models of the ionization effects. For a hot low-Z 
plasmas, such as a deuterium-tritium plasma during ICF ignition, the ions are likely to be fully ionized in 
any event, so this is not a serious restriction for clean thermonuclear burn with low-Z impurities. 
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employ a capital-I subscript to denote properties that correspond to the collective set of ions, 
such as the total ion number density n r = X]j n « or a common ion temperature T : . When 
equilibrium distributions are required in calculations, they are assumed to be Maxwell- 



Boltzmann, a 
more effort 



though a generalization to Fermi-Dirac statistics can be accomplished with 
3j. For problems involving hot thermonuclear burn, however, the fugacity is 
small and Maxwell-Boltzmann statistics is an accurate approximation. 

Let dS a b/dt denote rate at which the energy density of plasma species a changes because of 
its Coulomb interactions with another species b (the rate from the a-species to the 6-species). 
This rate is proportional to the temperature difference, and can be expressed by 

d -^ = -C ab (T a -T b ) . (f.l) 
at 

The sign convention in (11.11) implies that when the rate coefficients C ab are positive, then 
energy will flow from the hotter species to the cooler species, as it must. Section 12 of 
BPS used dimensional continuation to calculate the general rate coefficients C a b in a weakly 
coupled, but otherwise arbitrary, plasma. For simplicity, we will not perform the general 
calculation until Lecture III. In this and the following lecture, we will concentrate on the 
energy exchange between electrons and ions only. Since the electron mass m e is so much 
smaller than a typical ion mass m 1 , the electrons will come into equilibrium first with tem- 
perature T e on some time scale r e . The energy transfer rate among ions is a factor ^Jm e /m l 
slower than the corresponding rate for electrons, and therefore the ions will equilibrate 
to a common temperature T : in a time r : ~ r e ^Jm l /m e . Finally, as the electrons and 
ions exchange energy through Coulomb interactions, these systems too will equilibrate on a 
time scale r eI ~ r : ^m l /m e ~ r e {m l /m e ). Consequently, one finds a hierarchy of time scales 
r e -C Ti <C T e i, and it indeed makes sense to consider the electron and ion systems as having 
distinct temperatures T e and T : , with subsequent equilibration between them. Taking a = e 
and b = i in (jl.ip . and since the ions have a common temperature 7$ = T I5 the rate equation 
of interest is obtained by summing over the ion components of (11.11) to give 

^ = -C eI (T e -T : ) , (1.2) 
at 

where C eI = Ylii^ei and d£ el /dt = d£ ei /dt. 

The coefficient C el is the quantity we wish to calculate in this and the next lecture. This 
coefficient contains the energy-exchange physics between electrons and ions resulting from 
mutual Coulomb interactions, including possible collective effects and large-angle collisions. 
General expressions for the individual C e i were calculated in Section 12 of BPS. They are 
somewhat complicated and involve various one-dimensional integrals that can only be per- 
formed numerically. However, the collective rate coefficient C el simplifies considerably when 
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the mild restriction m e /T e <C m l /T 1 is imposed (a sum- rule is employed in the approxima- 
tion, and the simplification occurs only for C el = Yli^ei an d not for the individual C e i). If 
the high temperature limit is further imposed, then the rate can be written in a quite simple 
analytic form. 

For m e /T e <C m I /T l and T e l ^> e H : 

C eI = — uj x \hr-=r lnA BPS , with lnA BPS = - 

Z7T V Z7T 1 f, I 



, 3 8T " 



-7-1 



;i.3) 



where 7 = 0.57721 • • • is the Euler constant, n e and u e are the electron Debye wave number 
and plasma frequency, and u 2 = J2i u i * s sum °f the squares of the ion plasma frequencies. 
Since the small binding energy e H = 13.6 eV of the hydrogen atom sets the temperature scale, 
and since the condition m e /T e <C m I /T 1 is not very restrictive, the rate coefficient (jl.3p is 
applicable in almost all circumstances of interest. We shall devote the next lecture to deriving 
this expression. For now, note that equation (11. 3ft corresponds to Eqs. (3.61) and (12.12) 
of Ref. [l|, where I have taken this opportunity to correct a small transcription error: when 
passing from Eq. (12.43) to Eq. (12.44) in Ref. [l|, a factor of 1/2 was dropped. Restoring 
this factor of 1/2 changes the additive constant outside the logarithm from the —7 — 2 that 
appears in Eq. (12.12) of Ref. [l| to the constant —7 — 1 in (11.31) above. 

For reasons to be discussed shortly, rationalized units are preferred for dimensional con- 
tinuation, and I will employ this choice in all that follows. Nonetheless, expression (II .3p is 
written in a manner that does not depend upon this choice: the Debye wave number K e , and 
the plasma frequencies u e and uj^ can be calculated in your favorite units. For example, in 
Gaussian units where the electric potential takes the form V = e 2 /r, the Debye wave number 
and the plasma frequency of species b are given by = inelrib/Tb and u>l = 47r elnb/mb. 
In the rationalized units employed here, the electric potential is given by V = e 2 /Airr, and 
we have 

-Lb 

ut= e J^. (1.5) 

The square of the total Debye wave number is k 2 ^ = J2b an d t ne total Debye wave length 
is A D = k" 1 . 
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B. The Problem with the Problem 



Let us now consider an arbitrary plasma component of mass m, which I will otherwise 
leave unspecified, and let /(v, t) denote the Boltzmann distribution for this species. Then 
the average (kinetic) energy density of this component is 

£ = J d 3 v^mv 2 f(v,t) . (1.6) 

If we work to leading and next-to-leading order in the number density, calculating the energy 
exchange between plasma components will then involve keeping a tally only of the kinetic 
energy, as in (II. 6p . This is because the potential energy is higher order in the number 
density [or more precisely, the potential energy is higher order in the plasma coupling g, to 
be defined later in (|4.2p ] . As the system interacts with other plasma components through 
mutual Coulomb interactions, it will loose or gain energy depending on the temperature 
gradients with other species, and the energy exchange rate is given by 

§ = /'*' s™" < L7 > 

In contrast to (11. ip and (11.21) . for ease of notation I have temporarily dropped the plasma 
component subscripts on the rate, and I will keep with this convention until the final calcu- 
lation presented in Section IVIBI We see that the entire problem is bound up in calculating 
the rate of change df/dt from an appropriate kinetic equation that captures the relevant 
physics. As it turns out, however, there is a serious problem in performing all such calcula- 
tion with the Coulomb potential in three dimensions: the integrals in the kinetic equations 
diverge logarithmically, and they do so at both large and small distances. 4 

For processes in which large-angle scattering is important, such as the charged particle 
stopping power, it is natural to use the Boltzmann equation, which I will write in the 
abbreviated form 

f + V'V/ = B^ (1.8) 

where V is the gradient in position space, and B[f] is the scattering kernel, whose precise 
form will not concern us until the next lecture. The gradient vanishes because of spatial 
uniformity, so we will set v ■ V/ = 0. The Boltzmann equation was designed to account 
for the statistical effects of short- distance collisions, and although its original context was 
classical, quantum two-body scattering effects can easily be incorporated. In fact, since the 

4 It is curious that this problem occurs for the Coulomb potential only in three dimensions, the case of most 
physical relevance in plasma physics. 
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scattering phase shifts 5e are known analytically for the Coulomb potential, Ref. [l| used this 
to calculate the two-body quantum corrections to all orders [in the quantum parameter t] 
to be defined in (14. 141) ]. The kernel B[f] therefore contains all short-distance or ultraviolet 
physics, for both classical and quantum scattering. However, (in three spatial dimensions) 
the Coulomb potential is long-range, and the integrals in B[f] diverge in the infrared; or 
equivalently, if we write the scattering kernel B[f] in terms of momentum integrals, the 
divergence appears at small values of momentum. This was not a problem in Boltzmann's 
original formulation of (11.81) . since the Coulomb potential was unknown at that time, and he 
modeled particle collisions in terms of hard-sphere scattering. In summary, the Boltzmann 
equation gets the short-distance physics correct, including quantum two-body scattering, 
but it misses the infrared physics. The fact that the Boltzmann equation misses the long- 
distance physics manifests itself as an infrared divergence in the scattering kernel B[f], 
thereby rendering calculations meaningless (unless we tame, or regularize, this divergence). 

Given that the Boltzmann equation misses the long-distance or infrared (IR) physics, we 
might be tempted to try the Lenard-Balescu equation, which I will write in the abbreviated 
form 

?L+v.Vf = L\f\, (1.9) 

where L[f] is a scattering kernel whose exact form will be needed only in the next lecture. 
Again, the gradient term will be set to zero because of spatial uniformity. The Lenard-Balescu 
equation takes the form of a Fokker-Plank equation, with the kernel L[f] chosen to capture 
the correct IR physics. However, for the Coulomb potential (in three spatial dimensions), 
the Lenard-Balescu equation misses the short- distance or ultraviolet (UV) physics, and this 
is manifested by a UV divergence in L[f}. The situation for the Lenard-Balescu equation 
is exactly reversed compared to that of the Boltzmann equation. This is what Ref. [4J calls 
the "complementarity" of these two kinetic equations, and in the dimensional continuation 
procedure we will use this fact to our advantage. 



II. TRADITIONAL METHODS 
A. Heuristic Models 

The rate equation (11.71) reduces to a one-dimensional integral over the entire range of 
physical length scales (or momentum scales, if one so chooses), from zero all the way to 
infinity. Trouble arises for the Coulomb potential in three dimensions since the integral in 
question is logarithmically divergent at both integration limits. We must therefore regulate 
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the integral in some manner. Dimensional continuation is one such procedure, but there are 
others. This divergence problem was first worked around by simply cutting off the divergent 
integrals by hand, with the cutoffs themselves being chosen by physical arguments (rather 
than a calculation). 

The energy exchange rate we are considering is but an example of a larger class of problems 
involving characteristic, but disparate, length or energy scales in which the measured quan- 
tity of interest is (logarithmically) insensitive to the physics above and below these scales. 
For these problems, the simplest and most intuitive regularization scheme is to replace the 
offending integration limits, i.e. infinity and zero, by the finite and non-zero physical scales 
of the problem. These two scales then act as formal integration cutoffs, giving a finite loga- 
rithm of the ratio of the scales. Furthermore, because the system is insensitive to the physics 
above and below the respective cutoffs, this procedure provides a physically meaningful re- 
sult. Expressed in terms of length, we will denote the long- and short- distance scales by 6 max 
and 6 min , and the integral over scales leads to a finite logarithm involving the ratio of the 
physical length scales, so that 



where K is an easily determined prefactor with dimensions of energy density per unit time. 
As we shall see, a calculation to leading order in the number density is sufficient to provide 
the coefficient K, while a next-to-leading order calculation is required to find the exact terms 
under the logarithm. 

The problem with this regularization prescription, which I will call the heuristic scheme, 
is that we can only estimate the values of the physical scales b max and fe m j n to within factors 
of order one or so. For example, it is physically reasonable that the scale of the long 
distance cutoff in a plasma is set by a Debye length, so that fc max = ck~ 1 with c being 
a dimensionless constant of order unity; but what determines the exact value of c? In 
fact, how does 6 max ~ k^ 1 arise naturally from the kinetic equations themselves, rather than 
simply being chosen by hand? And should one use the total Debye wave number k u , or just 
the contribution from the electrons K e l The origin of the short distance cutoff 6 m i n is even less 
clear. In the extreme classical limit, we expect this scale to be set by the classical distance of 
closest approach r c \ between two colliding particles, so that 6 min ~ r c i. 5 On the other hand, we 
expect quantum effects to dominate when r c \ becomes smaller than some thermal Compton 
wavelength r qm ~ h/q, with q being a typical thermally averaged momentum transfer. In 

5 A further ambiguity arises in precisely denning r c \, as this involves a somewhat arbitrary choice in the 
thermal averaging procedure. 




(2.1) 
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this case we expect b m i n ~ r qm . Worse yet, the intermediate region in which neither classical 
nor quantum processes dominate is often realized in a weakly coupled plasma, and in this 
case it is even more obscure how one should choose 6 m ; n . We must interpolate between the 
extreme classical and the extreme quantum scales, but in the literature the exact procedure 
is always somewhat ad hoc. For example, a common choice is to define 6^ in to be the sum 
of the squares of the classical and quantum scales, so that 



I will have more to say about this in Section TlV A 31 

The heuristic scheme forces us to choose the specific forms of 6 max and b m \ n motivated 
by imprecise physical arguments or heuristic exercises, which leads us into the art of model 
building rather than systematic calculation. Indeed, the very notation that we must choose 
a cutoff is misleading, since the physics itself must conspire to render all integrals finite. 
Consequently, the heuristic method suffers from an unknown coefficient under the logarithm, 
and only the approximate value of the ratio 6 max /frmin can be determined with this method 
(in fact, this ratio varies across an order of magnitude over models in the literature, rather 
than factors of two or three). As we shall see, determining the constant under the logarithm 
exactly is equivalent to determining the next-to-leading order term exactly; therefore, models 
of the form (12.11) are accurate only to leading order, and no better. 

B. Convergent Kinetic Equations: Traditional Approach 

Rather than merely regulating the integrals in a rate derived from the kinetic equations, 
as with (12. ip . a more sophisticated approach involves regularizing the divergences in the 
kinetic equations themselves. In other words, the theory itself is regularized, rather than a 
particular quantity being calculated within the theory. The method of dimensional contin- 
uation falls into this category, albeit with somewhat more subtle mathematical machinery 
than traditional approaches. These approaches, of which Refs. are good examples, 

are summarized and placed into a common framework by Aono in Ref. 0]. As discussed in 
Section II Bl of this lecture, one can view the Boltzmann and Leonard-Balescu equations as 
providing complementary physics since they both succeed and fail in complementary regimes. 
The Boltzmann equation gets the short-distance physics correct, while the Leonard-Balescu 
equation captures the long-distance physics; conversely, Boltzmann and Leonard-Balescu 
miss the long- and short-distance physics, respectively. This complementarity motivates a 
class of kinetic equations of the form {^J 





dt 



+ v .Vf = B[f) + L[f)-R[f) , 



(2.3) 
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where R[f] is a carefully chosen "regulating kernel" designed to subtract the long-distance 
divergence of the Boltzmann equation and the short-distance divergence of the Lenard- 
Balescu equation. At the same time, the kernel R[f] must preserve the correct short-distance 
physics of the Boltzmann equation and the correct long-distance physics of the Lenard- 
Balescu equation (a minimal requirement of the regulating kernel R[f] is that it take a 
Hippocratic Oath to do no harm, at least to subleading order in the plasma coupling). 
Each term on the right-hand-side of (12.31) is separately divergent, but collectively they lead 
to a finite collision kernel if properly interpreted. 6 As an example, Gould and DeWitt 0] 
regulated the long-distance divergence of the Boltzmann equation by simply replacing the 
Coulomb potential with a Debye screened Yukawa-like potential, 



where, as not to confuse symbols, I write c as the base of the natural logarithm and e as the 
electric charge. While Ref. [?| performed this operation by hand using physical arguments, 
one could easily introduce a kernel R[f] to do the same. 

While the approach to convergent kinetic equations described by Aono might appear to 
be more rigorous than the aforementioned model building approach of (12.11) . it is no more 
systematic: methods based on (12. 3j) or its equivalent do not contain the ability to estimate 
their own error, i.e. they cannot determine their domain of applicability. There is nothing in 
the formalism of (I2.3P that keeps track of the plasma coupling constant, or the order to which 
we are working in this constant. Indeed, one does not generally think of (12.31) in terms of a 
perturbation theory. In contrast, Ref. [lj calculates the rate using a systematic expansion in 
the plasma number density, or more precisely, in a dimensionless plasma coupling parameter 
g [to be defined by (14.21) and discussed at length in Section HV A] , Although written in a 
disguised form, the BPS rate coefficient (11. 3ft is an expansion to leading and next-to-leading 
order in the coupling parameter: the leading order term goes like g 2 \ng, the next-to-leading 
order is proportional to g 2 , and the 0(g 3 ) term provides an estimate of the error. Translating 
the work of Gould and DeWitt [7] into the language of a perturbative expansion in a plasma 
coupling constant, it turns out that their result is valid to order g 2 and is in agreement with 

6 By "properly interpreted" I really mean that each term on the right-hand-side of (|2.3p should be separately 
regularized in some manner, rendering the individual kernels B[f], L[f], and R[f] separately finite: B[f] 
must be regulated at long-distances, L[f] at short-distances, and R[f] at both long- and short-distances. 
If this procedure is executed correctly, so that the long-distance divergences of R[f] and B[f] cancel, as do 
the short-distance divergences of R[f] and L[f], both in a consistent fashion, then the entire right-hand- 
side of (|2.3p remains finite as the cutoffs are removed. If R[f] does not disturb the 0(g 2 ) physics, then the 
convergent kinetic equation will be accurate to 0(g 2 ). I will have more to say about how one regulates 
long- and short-distances consistently in Section |Vj and how this can be a quite non-trivial process. 




(2-4) 



9 



BPS to this order. 7 However, in their final result, Gould and DeWitt retain spurious higher 
order terms in g 3 . I call these terms "spurious" because Ref. 0] did not calculate the full set 
of order g 3 terms, but only some of these terms. Indeed, the notion of a systematic expansion 
in a small dimensionless parameter does not enter their calculational framework. As this 
example shows, convergent kinetic equations can be more accurate than the heuristic model 
building technique of the previous section, but one cannot be sure of their accuracy until a 
comparison with a systematic calculation has been made, as with BPS and Ref. Q. 



III. BESSEL FUNCTION EXAMPLE 

A. Analogy with Dimensional Continuation 

I will illustrate the main points of dimensional continuation with an example involving 
the modified Bessel function Kq(x), with an emphasis on analytic continuation and how this 
can be used to extract leading and next-to-leading order behavior. This example was first 
presented in Ref. jsj], and for pedagogical purposes it was also included in Appendix A of 
Ref. [1]. This example contains all the essential features of dimensional continuation, but in 
a mathematically simple form, and while it is an imperfect analogy, as all analogies are, it 
is explicit in all its details. We will show that the modified Bessel function K (x) has the 
expansion 



K (x) =^-lnx j + ln2-7 + 0(x 2 ) = -ln( — x ) + 0(x 2 ) , (3.1) 




NLO 

to leading order (LO) and next-to-leading order (NLO) in x, with 7 = 0.577216 ■ ■ • denoting 
Euler's constant. This expansion is quite accurate for small values x, with an error of order 
x 2 rather than x for symmetry reasons. The asymptotic expansion (13.1 1) is a well known 
result but it is rather difficult to prove by conventional methods because of the non- 
analytic leading-log behavior. However, the method of dimensional continuation allows us 
to derive this result rather easily. The price one pays for this ease of derivation is that one 
must learn (or recall) a bit of mathematical machinery which, at first sight, seems unrelated 
to the problem at hand. 

We start with the general integral representation of the modified Bessel functions [1( 



7 For a detailed treatment of Gould and DeWitt in the context of BPS, see Appendix B of Ref. 
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As the notation in (13. 2p suggests, we can think of v as the dimension of space and the 
integration variable k as the wave number. In this analogy, the argument x corresponds to 
the dimensionless coupling parameter g of the plasma. The following dictionary provides a 
useful mnemonic in relating this mathematical example to the plasma physics problem of 
real interest: 

x — > plasma coupling g 
v — > spatial dimension 
k — > wave number 
dkk"^ 1 — > integration measure d u k . (3.3) 

Pushing our physical analogy further, if we think of v as being the dimension of space, 
then it should always be a positive integer, which I will express by the conventional set 
theory notation v G Z + ; however, nothing per sein the integral representation (13. 2p requires 
that v G Z + . We can therefore think of v in expression (13.21) as being a continuous real 
variable [y G E), or indeed, a complex variable [y G C) if circumstances warrant. Similarly, 
for a real physical system written in the appropriate integral form, there is nothing in any law 
of physics that prevents us from interpreting the dimension of space v as being a complex 
number. Continuing v from the positive integers into the complex plane is what I mean 
by "dimensional continuation." There will be times when we restrict our attention to the 
real numbers only, rather than the complex numbers in general, and I will refer to this as 
dimensional continuation as well. As we shall see, this procedure of taking v G M. or v G C 
will allow us to regulate otherwise infinite integrals in a systematic and perturbative fashion. 
Finite manipulations can then be performed, the divergent poles will cancel from physically 
measurable quantities, and afterward we can take v to the appropriate integer dimension (in 
this analogy we take v — > 0, rather than v — » 3 as we do for the physics problem). 

B. Leading Order Terms 

Let us first calculate the leading order in x behavior of K v (x) for positive and negative 
values of v. For small positive values of x, the leading order x-behavior can be obtained by 
replacing the exponential in (13 .2p by one, except in the regions k — > and k — > oo, where the 
exponential is required for convergence. Taking v < first, note that the integral (13. 2j) is 
dominated by small values of k near the lower limit of integration. In terms of the analogy 
(13. 3p where k is a wave number, this corresponds to the situation in which long-distance IR 
physics is dominant. Therefore, when v < and the integral is dominated by small values of 



11 



k, the leading order contribution to (13. 2p can be obtained from the leading order behavior 
of the exponential, that is to say, the replacement 

exp|-| (k + ~jj expj-^-j with |x| < 1, z/<0 (3.4) 

will capture the entire leading order in x behavior for negative values of v. Note that 
(13 .4p provides convergence as k — > 0, while large- convergence is provided by the prefactor 
k v ~ x since v < 0. We will denote this leading order contribution by K<(x), and using the 
substitution (13.41) we write 

K<(x) = - / dkk v ~ l e- x/2k = ~ (~ ) T{-v) . (3.5) 



2 J 2 V2 

In the last equality of (13.51) . we have made the variable change y = x/2k, and we have used 
the standard integral representation for the Gamma function, 

POO 

T(z)= / dyy z ~ l e-y . (3.6) 



This representation of the Gamma function only converges when z > (actually, (13. 6p is 
valid over the complex z-plane with Rez > 0). When z = — v and v < 0, we may indeed use 

TO- 

We can find the leading order in x contribution when v > in a similar manner. In this 
case, the integral is dominated by large values of k at the upper limit of integration, and we 
make the substitution 

exp \~ ( k + ~ J 1 -> exp \-%- \ with |ar| < 1 , v > , (3.7) 



2 



thereby giving the leading order result 



1 r°° ,1 /r\-v 

K>{x) = - / dkk^e-^ = - Q T(y) . (3.8) 

Note that the exponential provides convergence as k — > oo, while the integrand possesses an 
integrable singularity at k = for < v < 1 (the integrand is non-singular at k = when 
i/ > 1). 

We will eventually take the v — > limits of (13.51) and (13. 8p . since we are interested in 
Kq(x) and not K u (x). While we could take v to be a general real or complex number 
until the limit is required, it is easier to consider only small values of v from the start (we 
need work no higher than linear order, since this and higher orders vanish when v — > 0). 
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To do this, we expand the Gamma function T(±z/) to linear-order in its argument using 
T(z) = 1/z — 7 + O(z). Taking z = ±v in this expansion gives 

id) 
= "i © 

Expressions (13. 9p and (13.101) are accurate to linear order in the dimension u; on the other 
hand, (13.91) gives the leading order in x contribution to K v (x) as defined by (13. 2\ when u > 0, 
and (13.101) gives the leading order in x contribution to K v (x) when v < 0. To compare these 
two expressions, we must analytically continue them to a common dimension. We will discuss 
this further in the next section. 

In terms of our physics analogy, expression (13.91) captures the leading order short-distance 
physics in the v > regime; the pole at v = corresponds to a small-fc divergence, which, 
pushing our physical analogy again, would reflect missing or incomplete long-distance physics 
(as with the Boltzmann equation). The situation is completely reversed for (I3.10p . which cap- 
tures the leading order long-distance physics for v > 0, with the pole at v — corresponding 
to a large-fc divergence arising from missing short-distance physics (like the Lenard-Balescu 
equation). As functions of z/, we see from (13.91) and (I3.10p that K>{x) and K<{x) are an- 
alytic in u, except for the simple pole at v = 0. As we shall see, the analytic continuation 
to complex v takes the same functional form as the individual expressions (13.91) and (I3.10p . 
each defined separately for v > and v < 0, respectively. 



1 — b"~f 



LO in x when v > 



(3-9) 



LO in x when v < 



(3.10) 



C. Some Comments on Analytic Continuation 

Since analytic continuation plays such a central role in dimensional continuation, at least 
mathematically, I would like to briefly discuss the conditions under which a function can be 
analytically continued from one region of the complex plane to another. Recall that a function 
/ is said to be analytic at a point zq in the complex plane C, if and only if its derivative exists 
not only at zq, but also at every point within some open neighborhood of Zq. A function 
/ is analytic on a domain D in the complex plane if it is analytic at each point in D. 
Analyticity is a very stringent condition on a function, since the existence of the derivative 
of a complex function is a much more robust constraint than the corresponding existence 
of the derivative of a function on a real domain. This is because in the two-dimensional 
complex plane, the limiting procedure defining the derivative must exist regardless of the 
direction used in taking the limit. In fact, analyticity at a point zq is so strong that it implies 



the existence and continuity of all derivatives f^ n \zo) for any order n > ll|. In other 
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words, an analytic function on D can be thought of as being infinitely smooth on D, even 
though the definition of analyticity itself invokes only the existence of the first derivative, 
albeit on a neighborhood. 

Analyticity is such a stringent condition, that the behavior of an analytic function in a 
small domain is enough to determine its behavior in a larger region. Even if the function 
is only known along a one-dimensional curve in the complex plane, such as a portion of the 
real axis, this is enough to uniquely determine the function in the complex plane. 8 This is 
because the derivatives of an analytic function exist to all orders, and therefore a Taylor 
series expansion about a point of analyticity always exists with some non-zero radius of 
convergence. For example, one of the powers of analytic continuation is that the original 
function may be defined in any manner over a subregion, and this can generate a unique 
function over a larger region. As an application of this, take the analytic function defined 
by the infinite geometric series 

oo 

/(z) = X>, (3.11) 

n=0 

a series that converges only for \z\ < 1. Upon defining / by (13.111) . we therefore take the 
domain D x to be the unit disk about the origin, excluding the unit circle itself. In this 
domain, the geometric series converges to 

oo - 

* n = y— ^ f o r zeB 1 . (3.12) 

n=0 

Note, however, that the function g(z) = 1/(1 — z) is defined over the entire complex plane 
except z — 1, a region I will call Ci. Since the function / is only defined within the unit 
circle, and since / and g agree within the unit circle, the function g:C\ —>■ C is the unique 
analytic continuation of / : D\ — > C. 

As a more relevant example, consider K<{x) with v < as given by (13.101) . To compare 
this with K>(x), which is determined by (13. 9p for v > 0, we must analytically continue 
K<(x) to the positive real u-axis. We can think of K<{x) as a sequence of functions of 
an independent variable v indexed by a continuous label x; therefore, in a more suggestive 
notation, we temporarily write f x {y) = K<(x) . While the collection of functions f x {v) are 
defined in (13.101) on the negative z/-axis 1R~ (excluding the simple pole at zero), they can be 

8 In fact, Carlson's Theorem [l^] can be used to uniquely extend a function defined only on the integers to 
the whole complex plane. This is actually the theorem of most relevance here; however, as we shall see 
in Section IIVB1 in practice we can analytically continue a function defined on the integers without this 
Theorem. We appeal to Carlson's Theorem only to guarantee the uniqueness of this procedure. 
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analytically continued to the complex z/-plane C. Furthermore, the functions take the same 
algebraic form on the complex plane, namely, 



Uv) = ~h (I 



1 + ^7 



for v G C . 



(3.13) 



We can now restrict our attention from C in general to the positive z/-axis 1R + (excluding 
zero). This allows us to directly compare K<{x) and K>(x) at v G M + using the same 
algebraic forms as given by f 1 3 . 1 j) and (13.91) . In the next section, we discuss the implications 
of analytically continuing K<(x) from the negative axis v < to the positive axis v > 0. 
Alternately, we could continue g x {y) = K>{x) to the region v G M~ using the same functional 
form as ( 13.91) . and compare this with K<(x). 



D. Next-to-Leading Order Term 



We now illustrate the key mathematical result that allows dimensional continuation to 
extract not only the leading, but the next-to-leading order terms. Recall from (13. 9p and (13.101) 
that K>(x) and K<(x) are both leading order in x for v > and v < 0, respectively. Since 
these functions were calculated for mutually exclusive values of u, one might think that they 
cannot be compared. When viewed as an analytic function in the complex z/-plane, however, 
we have seen that K<(x) is also a function over the domain v > 0, in which case both K<(x) 
and K>(x) can be compared at the same values of v and x. Since the algebraic form is so 
simple, K< (x) takes the same functional form when analytically continued to v > as it did 
for v < 0. As illustrated in Fig. [T], this means that K<(x) becomes next-to-leading order in 
x along the positive real axis: 



K>(x) 



K<(x) 



1 

27 
1 

27 



z/7 



1 + 2^7 



LO in x when v > 



NLO in x when v > 



(3.14) 
(3.15) 



The expressions (13.141) and (13.151) remain accurate to linear-order in v. 

To see that (I3.15P is indeed next-to- leading order in x relative to (I3.14p . note that the 
x-behavior of the leading order contribution for v > can be written K>(x) ~ aH^. I have 
used the absolute value \u\ to emphasize that the power of x in (13.141) is strictly negative 
when v > 0. Similarly, along the positive v-axis we find the behavior K<(x) ~ x^ for v > 0, 
and we see that x~^ for < x <C 1. This means K>(x) K<(x) for > and 

< x < 1, and we are therefore justified in calling K>{x) leading order in x and K<{x) 
next-to-leading order. 
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analytically continue 
around the v = pole 




K<(x) ^ K<(x) 

x -\u\ x \u\ 

LO: large when x <C 1 NLO: small when x<l 

FIG. 1: The analytic continuation of K<(x) from v € K~ to ^ € K + in the complex i/-plane: the same 
expression can be used for K<(x) throughout the complex plane since the pole at v — can easily be 
avoided, as indicated in the figure. Note that K<(x) ~ aH"! is leading order in x for v < 0. However, 
upon analytically continuing K<{x) to v > 0, the x-dependence becomes K<{x) ~ a;'"', and the function is 
next-to-leading order relative to K>(x) ~ a; - . 

Strictly speaking, we have only shown that K<(x) is subleading relative to K>{x) when 
f > 0. To conclude that K<(x) ~ a; 1 ' is indeed next-to-leading order relative to K>(x) ~ a; - ' 7 , 
it is important to establish that there are no powers of x between x v and x~ v in the expansion 
of K v (x). For v > 0, one simply subtracts (13.91) from (13.21) . and it becomes clear that this 
error is higher order in x than x v '. For a more detailed proof of this, see footnote 2 of 
Ref 



■ k 



A similar statement holds for K>(x), namely, as we analytically continue from 
v > to v < 0, the quantity K>(x) switches from leading order to next-to-leading order 
in x relative to K<(x). 



E. Assembling the Pieces 

We have now assembled enough results to find K (x) to leading and next-to-leading order 
in x: we simply add the expressions ( 13. 14ft and ( 13. 15ft and take the limit of vanishing v. Note 
that this does not lead to any form of "double counting." Instead, we are simply adding 
the next-to-leading order term (13.151) to the leading order term (13. 14j) at a common value 
of v > 0. Upon taking the limit of vanishing u, or more precisely v —>■ + since v is always 
positive in ( 13.141) and ( 13.151) . we obtain K (x) to leading and next-to- leading order in x. 

We now calculate this limit, proving that 



lim 

2^0 + 



K>(x)+K<(x) 



-lnx + ln2-7. (3.16) 



Let us first expand (x/2) ±u in powers of v. We will denote z = x/2, from which we find 
z ±u = e ln ( z± "^ = e ±ulnz = 1 ±vhiz + 0(u 2 ), or in summary: 



z ±v = l±v\nz + 0{v 2 ) . (3.17) 
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This allows us to express (I3.14p and (13.151) as 



and 



K<(x) 



1 

2^ 



i+ "Kf) 



l — 1/7 



1 + ^7 



+ 0(1/) 



+ 0(z/) 



2f 2 V2 



J--lln^ 
2u 2 V2 



I + 0{v) , (3.18) 



7 



+ C?(i/). (3.19) 



When we divide (13. 17ft by a factor of u, as required by (13.141) and ( 13. 15ft . note that we find: 
(i) a pole 1/v from the first term in (13. 17j) . (ii) a non-analytic finite contribution ±lnz from 
the second term, and (iii) the error in v becomes 0(v), which is the same order in v that 
we are neglecting in (13. 14j) and (I3.15p . The error in u, which vanishes in the limit v — > + , 
should not be confused with the error in x, the latter being 0(x 2 ) for vanishing v. Note that 
the pole terms cancel upon adding (13 . 1 8[) and (13.191) . so that 



K>(x)+K<(x) 



ln(|)- 7 + OM 



(3.20) 



thereby giving (13.161) as v — > + . As we have discussed, there are no x-dependent terms that 
lie between K>(x) ~ x~ v and K<(x) ~ x u , so this procedure has captured the leading and 
next-to-leading order behavior in x. 

In exactly the same way, we can also calculate the leading order and next-to-leading order 
contribution to K (x) by taking the limit from the left, 



lim 



K<(x)+K>(x) 



In x + In 2 — 7 



(3.21) 



I should point out a potential notational problem in (13.161) and (13.211) . Concentrating on 
(I3.16P for the moment, the limit v — > + indicates that both terms K> and K< are understood 
to live in dimensions v > 0, with the second term K< having been analytically continued 
from v < 0. The notation with which the term K< is written in (13.161) . however, does not 
indicate that it has been analytically continued. This should be no cause for confusion, 
however, since K< takes the same functional form in any dimension v\ therefore, a separate 
notation indicating that the K< in (13. 16[) has been analytically continued is unnecessary. We 
can simply add K> and K< as calculated in v > and v < respectively. 
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IV. DIMENSIONAL CONTINUATION 



A. Rate of Energy Exchange as a Perturbative Expansion 



1. The g- Expansion 

Before moving on to the details of dimensional continuation, we must first discuss the 
plasma expansion parameter g. Since the problems in plasma physics are usually so compli- 
cated as to preclude a perturbative approach, most plasma physicists do not usually think in 
terms of expanding systematically in a small dimensionless parameter. However, for a weakly 
to moderately coupled plasma, it is a quite fruitful approach to perturbatively expand the 
rate in a small dimensionless coupling constant. 



That such a universal parameter for a plasma exists was discussed at length in Ref. [131 ]. 
where it was shown that any physical quantity associated with a plasma whose species are 
in equilibrium with themselves (such as the plasma we are studying) can be expanded in 
integer powers of a dimensionless coupling constant g defined by 

_ Coulomb Energy for Two Charges at Separation A D (4 1) 

Temperature in Energy Units 

with A D being the Debye length of the plasma. Since the potential energy between two like 
charges is given by V — e 2 /4nr in rationalized units, and writing the Debye wave number 
clS = Aj, 1 , the coupling parameter is therefore 

' = £f ■ (4 - 2) 

For a multicomponent plasma, there is actually a coupling constant for each pair of compo- 
nents, 

g ab = -: „ , 4.3 
ATcT b 

with Kb defined by (11.41) . However, when the pairs have approximately the same coupling 
strength, then the single parameter (14.21) adequately characterizes the entire plasma. Ex- 
pressing the charges as e a = Z a e and e& = Zb e, we can write g a b as 

1 b 

and we see that the coupling constant is proportional to the cube of the electric charge, the 
square root of the density, and the inverse (3/2)-power of the temperature. 
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Recall that the usual plasma parameter T is defined in a similar manner to (14. ip . except 
that the charge separation is determined not by A D , but by the inter-particle spacing in the 
plasma, 

j_, _ Coulomb Energy for Two Charges at Inter-particle Separation d ^ 

Temperature in Energy Units 

The inter-particle spacing d is defined in several ways in the literature, but the idea is to 
transform the plasma number density n into a length scale, so that d oc n -1 / 3 . The most 
common convention is to define d to be the radius of a sphere containing, on average, a single 
plasma particle, so that And 3 / 3 = 1/n, and therefore 



e 2 AtttiN 1/3 



4vrT 



(4.6) 



With this convention, the relation between the two plasma coupling parameters for a single 
plasma species is T 3 = g 2 /3, and for an arbitrary number of plasma species we always find 
g oc T 3 / 2 . We can therefore use either g or T to characterize the strength of the plasma, as g 
and T become large or small together. 9 For our purposes, however, there is an an advantage 
to using g rather than T since physical quantities are expanded as integer powers of g, while 
they expand in fractional powers of T. 



2. Next-to- Leading Order and the Coulomb Logarithm 

As I have said, any plasma quantity can be written as a power series expansion in integer 
power of g, with the possible exception of non-analytic terms involving In g. For the process of 
energy exchange via Coulomb interactions, this non-analyticity arises from the competition 
between disparate physical length scales. As an expansion in g, the rate of energy exchange 
takes the form 

^ = - A^]ng + Bg^ + 0{g 3 ) . (4.7) 

LO NLO 

9 The dimensional continuation method requires that we perform calculations in an arbitrary spatial di- 
mension v; however, the parameter g — e 2 K u /4irT as given by (|4.2j) is dimensionless (in the engineering 
sense) only in three spatial dimensions. As such, it is meaningless to call g large or small in any dimension 
other than j/ = 3. This is no cause for alarm, however, since in footnote 12 11 I will construct a dimensionless 
coupling parameter g v in arbitrary spatial dimensions v. The parameter g v can then be used to charac- 
terize the plasma strength in v dimensions, and it will have the properties that (i) g v oc g and (ii) g v — » g 
as v — » 3. Property (i) implies that a (^-expansion is the same as a g-expansion, and (ii) implies this 
correspondence continues down to three dimensions, where g becomes the relevant expansion parameter. 
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In (14.71) . I have indicated the leading order in g (LO) and the next-to- leading order in 
g (NLO) terms in the (/-expansion: the first term is leading order relative to the second 
because |g 2 lng| > g 2 for small g. The minus sign on the leading order term of (14. 7p is a 
matter of convention. Since the logarithm Ing will be negative in a weakly coupled plasma 
(recall g < 1), the minus sign renders the coefficient A positive when the energy exchange 
is positive. The coefficient A was first calculated by Spitzer. The coefficient B, however, 
is very difficult to calculate, and this was the main purpose of BPS It is convenient to 
define the dimensionless coefficient C by B = —A In C, in which case we can write 

— = Ag 2 lnA coul + 0{g 3 ) , with lnA coul = -ln{Cg} . (4.8) 

We see, then, that knowing the next-to-leading order term is equivalent to knowing the exact 
coefficient C under the logarithm. Note that the minus sign renders the Coulomb logarithm 
positive when g is very small, in keeping with convention. 



3. Factors of g Inside the Coulomb Logarithm 

For the heuristic model building of Section III Al let us pause for a moment and show that 
the argument of the Coulomb logarithm in (12. ip is indeed proportional to g, as required by 
(14. 8p . On physical grounds we saw that the long-distance scale 6 max is set by a Debye length, 
and therefore we nominally set 6 max = k" 1 . In the extreme classical limit, the short-distance 
cutoff 6 min is set by the classical distance of closest approach r cl , so that b min = cr d . For 
simplicity, we will choose the coefficient c such that K = Ag 2 between (12. ip and (14.81) . in 
which case 

lnA coul = ln(m . (4.9) 



Let us consider two unit charges of mass m approaching one another with zero impact 
parameter. The rms speed of each particle is determined by 



\mv-=\ T =► v=f-^, (4.10) 

while energy conservation \ mv 2 + \ mv 2 = e 2 /47ir c i gives the distance of closest approach, 

e 2 e 2 I 

rd = -: ~ 2 = y- ^ ■ (4.11) 

A7TmV z 47T 61 

In the extreme classical regime, we see that the argument of the Coulomb logarithm in (14.91) 
is indeed proportional to the plasma coupling constant, 

Cr d K D = o = ~9 ■ (4-12) 



b mxK 3 iTi T 3 
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Let us now look at the ad hoc interpolation (I2.2p between the classical and quantum 
regimes. Up to this point I have said very little about quantum mechanics. While I will 
not dwell on quantum corrections, I will briefly discuss a dimensionless expansion parameter 
that characterizes the strength of the quantum two-body scattering correction. There are 



many ways of defining such a parameter, but I will follow Ref. 13|, taking 



_ Classical Distance of Closest Approach 
Thermal Wavelength 

With this definition, quantum corrections are large when rj <C 1. Motivated by the de Broglie 
wavelength of a particle, the thermal wavelength of a plasma species is given by r qm = h/q, 
where q = mv is a typical momentum transfer suffered during a collision. Definition (14.131) 
yields rj = r d /r qm = (e 2 / 4nmv 2 ) ■ (mv/h), or more succinctly 

AttHv 



V = J^-=.i (4-14) 



from which (12.21) gives 



1/2 

1 + - ) g. (4.15) 



3 V V 



In the extreme quantum limit in which rj -C 1, this becomes 6 m in/Vax = (c/3) (g/rj). 

Finally, note that the factors inside the BPS Coulomb logarithm (11. 3B are also proportional 
to the coupling constant g, and upon dropping the electron subscripts for convenience, we 
find 

(4.16) 

T 2 3 rf K ' 

We see that the ^-dependence of the Coulomb logarithm arises quite naturally. However, 
the accompanying coefficient under the logarithm might also possess ^-dependence, thereby 
obscuring the (^-dependence unless we are careful. 



B. Mathematics of Dimensional Continuation 

Before describing what dimensional continuation is, allow me to first state what it is not. 
Dimensional continuation is not performing an integral to a fractional power of the spatial 
dimension, as with the meaningless expression 

oo 

dhf(k) . 

o 
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Instead, dimensional continuation is the following. Suppose some physical quantity of interest 
can be written as an integral over a kernel 



where / is determined by the physical equations of motion, whether classical or quantum. 
The integrand is of course a function of the physical parameters, such as the masses and 
charges of the fundamental particles, and I have abbreviated this dependence by the param- 
eter m. For definiteness, we will think of A; as a wave number with dimensions of an inverse 
length. The laws of physics, from which (14.171) follows, are usually written in three dimen- 
sional space. Thus, we usually take k to be a three-dimensional vector, and we integrate k 
over the entire three-dimensional Euclidean space M 3 . 

The known fundamental laws of physics themselves, however, do not specify a particu- 
lar spatial dimension in which they hold. In fact, as far as the known laws of physics are 
concerned, the actual value v of the spatial dimension can be viewed as a free integer param- 
eter: it is simply an unexplained empirical fact that we live in three dimensions. 10 We can 
therefore express any three-dimensional physical quantity or law, such as Newton's equation 
of motion or Gauss' law, in any number of integer dimensions v G Z + . For example, we can 
write down a corresponding quantity to (I4.17P in an arbitrary number of dimensions, 



where the wave vector k is now a //-dimensional vector, and we integrate over the entire 
//-dimensional Euclidean space W . I have placed a subscript on the integrand f v to indicate 
that it is determined by the theory expressed in v dimensions. At this point, the spatial 
dimension is a non-negative integer, so that v = 1, 2, 3, 4, • • ■ . Since the integer v is arbi- 
trary in (14. 181) . I have indicated that the corresponding quantity contains //-dependence by 
writing Q(u; m); however, for notational simplicity I will drop the parametric dependence of 
quantities such as mass and simply write Q{v)- A quantity that diverges in three-dimensions 
will be finite in arbitrary z/, but it will typically exhibit a simple pole of the form 



One would expect that the fundamental theory of nature, a theory of everything, would predict the number 
of space-time dimensions, solving this mystery at last. One of the great successes of string theory is that it 
is one of only two known theories that indeed predicts the number of space-time dimensions - the theory 
is inconsistent in all but nine space and one time dimensions. Accordingly, this is also one of the great 
failures of string theory. The other theory is super-gravity, which is only consistent in eleven dimensions. 




(4.17) 






(4.19) 
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where Qi(v) is finite at v = 3. 11 Dimensional continuation is simply the act of treating Q(v) 
as a function of a complex argument u, after the integral (14. 18j) has been performed for all 
positive integer values of v. This is not an unfamiliar procedure, as the factorial function v\ 
with v G Z + can be generalized to the gamma function T(y + 1) in which v G C Indeed, 
given any function Q[y) defined on the integers, with only mild restrictions placed on the 
function at large values of the argument, Carlson's Theorem [l^] allows us to continue this 
function uniquely to the complex plane. 

How does dimensional continuation work in practice? We will look at a few specific 
examples in Sections IIVCI and IVTl but for now let us consider a general physical quantity Q 
in which the integrand in (14.171) depends solely upon the modulus of k, so that 

Q = J d 3 kf(k) . (4.20) 

In such a case it is not uncommon that the integrand in the generalization (14.181) is only a 
function of the modulus of the //-dimensional wave vector k, with the same functional form 
as the integrand in (14.201) . In other words, in (14.181) we have /^(k) = f(k) with k = |k|, 
thereby allowing us to write 



Q(u) = jd v kf(k). (4.21) 
Since the integrand is a function only of fc, we can extract the angular integrals and write 

n poo 

Q(v)= dto v -i ■ / k v ~ x dkf{k) . (4.22) 



o 



(y—1) integrals one-dimensional integral 

At this point, the dimension v is simply an arbitrary positive integer, v G Z + . As we will 
show in the next paragraph, the integration over all angles gives 

r 2W 2 

<Vi = J = --- . (4.23) 

This leaves a one-dimensional integral to perform, in which v simply acts as a parameter, 

F{v) = / dkk v ~ x f{k) . (4.24) 
Jo 

The physical quantity now becomes the product of (14.231) and (I4.24p with v G Z + . In the case 
of (14.231) . we already know how to analytically continue v to complex values. On the other 



11 As we shall see in the next section, it is the nature of the Coulomb force that gives a pole in three 
dimensions, rather than some other value of the dimension. We shall further see that there is important 
physics in the residue Qo of the pole. 
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hand, for (14.241) we can think of v as being an arbitrary complex number when performing 
the one-dimensional integral over k, and therefore we can regard F{u) as a function over the 
complex z/-plane, thereby giving 

Qiy) = F{u) with ueC. (4.25) 

In this manner, we can regard Q{v) as a function of a complex argument u, and by Carlson's 
Theorem 12], this is the unique continuation from positive integer values of v to complex 
values of v. 

As an example of this procedure, let us prove (I4.23p . First, consider the one- dimensional 
Gaussian integral 



CO 



dke~ k =v^F. (4.26) 
If we multiply both sides together v times (with v £ Z + ), we find 

/oo poo poo p 

dk ie - k i dk 2 e-^--- dk y e~ kl = d u ke- k2 , (4.27) 
-oo J-oo J —oo J 

where the wave vector k in the exponential of the last expression is the v- dimensional vector 
k = (hi, k 2 , ■ ■ ■ , k v \ and k 2 is the i/-dimensional inner product k 2 = X^=i • As in f |4.22j) . 
we can factor the angular integrals out of the right-hand-side of (14 .2 7ft . and the remaining 
one-dimensional integral can be converted to a Gamma function with the change of variables 
t = k 2 : 



71 ' 



//•OO p -I poo p -1 

dn^-J dkk u - 1 e- k ' 2 ''=jdn v -i- - J dtt vl2 - x e~ l =JdVL v _ l ■ -I>/2). (4.28) 



Solving for J dVt u ^\ in (14.281) gives (14.231) . As an aside, it is interesting to note that we have 
also found the hyper-area of a (u— l)-dimensional sphere of radius r in M. v , 

2 W 2 

K = j^j , (4.29) 
and integrating (14.291) gives the hyper- volume of a v- dimensional ball of radius r, 
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C. Physics of Dimensional Continuation 

So far we have only introduced dimensional continuation as a regularization prescription, 
a means by which infinite theories can be rendered temporarily finite. This is the use to 
which dimensional continuation is commonly employed in field theory; however, this alone 
is not sufficient to render the method suitable for our purposes. Instead, there are also 
physics reasons that make dimensional continuation applicable to the problem at hand. I 
will now show that the Coulomb potential, and indeed physics in general, behaves differently 
in different spatial dimensions. In particular, short-distance or ultraviolet (UV) physics 
dominates when v > 3; while long-distance or infrared (IR) physics dominates when v < 3. 
In v — 3 dimensions, both UV and IR physics are equally important. 12 This means that 
dimensional continuation acts as a "physics sieve," allowing us to capture the leading UV 
and IR physics simply by performing the relevant integrals in dimensions greater than or 
less than the traditional v = 3. In the next few paragraphs, we will discuss why changing 
the dimension of space emphasizes either long- or short-distant physics. Understanding this 
point is crucial for all that follows. 

We now turn to finding the //-dimensional Coulomb potential, an exercise that succinctly 
illustrates how the physics of a system changes with the dimension of space. Besides illus- 
trating that UV physics dominates in higher dimensions (and conversely), this example will 
also clarify the manner by which one performs physical calculations in arbitrary dimensions. 
Let us consider Poisson's equation 



There is nothing in this equation per se that restricts us to three dimensions. We choose 
to describe the electric field E and the spatial coordinate x as three-dimensional vectors be- 
cause we live in three dimensions. However, from the mathematics alone, we could equally 
well consider these vectors as living in an arbitrary //-dimensional space W. The rectilin- 
ear coordinates would then become x = (xi,X2, ■ ■ ■ ,x„), with a similar expression for the 
//-dimensional electric field, while the gradient would be V = (d/dx-y, d/dx 2 , • • • ,d/dx v ). 
It is convenient to write Poisson's equation in its integral representation, which will allow us 
to calculate the electric field of a point charge with relative ease. Consider a point-charge at 

12 There are a number of such coincidences, both physical and mathematical, that suggest there is something 
special about living in three spatial dimensions and one time dimension. This could well be the Anthropic 
Principle at work. 

13 One can most easily write the complete set of Maxwell's equation in a general dimension by employing the 
Lorentz covariant form ^2a=o®F a ^ /dx a — where the electric and magnetic fields have been expressed 
in terms of the anti-symmetric field tensor F a/3 . 




(4.31) 
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the origin given by p(x) = e<5 M (x), and let £ be any volume containing the charge e. Then, 
in a general number of dimensions, we can integrate (I4.3ip to obtain 

cfxV-E = e. (4.32) 



s 



Note that the dimensionality of space is now explicitly indicated by the integration measure. 

For our purposes, the advantage of the integral representation of Poisson's equation is 
that the electric field of a point charge in the v- dimensional space can easily be calculated 
by the same symmetry principles that hold in three dimensions. Suppose the volume X is 
a spherical ball B r of radius r centered on the charge. The boundary of B r is a sphere of 
dimension v— 1 and will be denoted by dB r (I am using the common notation d in differential 
geometry for the boundary of a manifold). From rotational symmetry, the electric field E 
of a point charge is directed radially outward and lies normal to the surface dB r at each 
point. We will denote the magnitude of the electric field at radius r by E(r). Recall that 
in (I4.29p . we calculated the hyper-area of the [y — l)-sphere dB r to be A u = Q u -ir u ~ 1 
with = u v l 2 /T(v/2). Since the divergence theorem holds in an arbitrary Euclidean 

space (like the laws of physics, there is nothing in the divergence theorem to restrict the 
dimensionality of space to three), we can write 

e = I ftV-E= I dA-E = Q u _ 1 r u - 1 ■ E(r) . (4.33) 
At position x, the electric field therefore takes the form 



where I am using the notation x = r x, with x being a unit vector pointing in the direction 
of x and r = |x| being the magnitude. It is often more convenient to work with the electric 
potential, a scalar quantity 4> v {r) defined by E{r) = —d(fi u (r)/dr. In fact, we need the 
potential energy V u (r) = e(f> u (r), so that 

Mr) = ^r£z, (4-35) 

with 

a = !),_ 1 (-2) = f ^ Iy . (4.36) 

For v = 3 we have C3 = 1/47T, which is the origin of the 4ir of rationalized units. Note 
from (14.351) that the engineering unit of electric charge is a function of the dimension v. 
This is because the i/-dimensional Coulomb potential (14.351) must have units of energy, and 
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t T 

UV r IR 

FIG. 2: Short-distance or ultraviolet (UV) physics dominates in dimensions v > 3. Long-distance or infrared 
(IR) physics dominates when v < 3. UV and IR physics are equally important in v = 3. 

consequently the engineering unit of e 2 is Energy x (Length)^ -2 = Massx (Length) l/ /(Time) 2 . 
It is quite natural that a composite quantity, made from the fundamental units of Mass, 
Length and Time, change its engineering dimension with the dimension of space. 

Figure [2] shows the Coulomb potential for z/ = 3, along with two representative dimensions 
on either side of v = 3. As the figure illustrates, the short-distance or UV behavior of 
the Coulomb potential becomes more severe as the dimension increases above v = 3, while 
the long-distance or IR behavior dominates for dimensions below v = 3. The arbitrary 
integration constant for the potential energy has been adjusted in each case so that all three 
graphs intersect at a single point. This was purely for aesthetics, as it renders the differences 
between the potentials more apparent. Despite the trouble we went through in the previous 
paragraph to find the coefficients C u , in this paragraph (and only in this paragraph) I 
have temporarily set C v — 1. This will make it easier to compare the r-dependence of 
various potentials, and I would rather opt for clarity over notational consistency. For the 
representative potential with dimension below v = 3, I chose to graph the one- dimensional 
potential V\(r) = e 2 r rather than the two-dimensional potential V^(r) = e 2 ln(r/ro), where ro 
is an arbitrary integration constant with units of length. In both cases the potential increases 
without bound at long distances, thereby illustrating the dominance of long-distance or IR 
physics in dimensions less than three. However, I chose to graph V\{r) rather than V^r) 
because the latter possesses a potentially misleading divergence as r — > 0: unlike the short- 
distance or UV divergence associated with the potential in dimensions greater than three, the 
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projectile v p t 




66 



FIG. 3: A projectile of charge e p , mass m p , and velocity v p passing a fixed charge e^. The impact pa- 
rameter b is normal to the velocity, so that b • v p = 0. The radial separation between the charges is 
r(t) = (b 2 + Vpt 2 ) 1 / 2 , and only the 6-component E ± — Ecosd — (b e^/O^-i) r~ v of the electric field con- 
tributes the impulse integration (|4.40|) . where cos6> = b/r. 



r — > divergence of V^r) is integrable, and consequently harmless. Therefore, for purposes 
of illustration, the potential Vi(r) = e 2 r makes the point better than V^(r) = e 2 ln(r/ro). 
Remarkably, we now see that by simply selecting the dimension u, we can dial a potential 
V v (r) that filters either long-distance or short-distance physics. 

I would now like to show how the pole at v = 3 arises from the v- dimensional Coulomb 
potential. In the temperature equilibration process, individual plasma species exchange 
energy by mutual Coulomb interactions. For example, consider a particle in the plasma with 
charge e p , mass m p , and velocity v p , and suppose it passes another charge with an impact 
parameter b. This is illustrated in Fig. [31 where, to zeroth order, the projectile follows the 
straight line x(t) = b + v p t as a function of the time t, with b • v p = 0. 14 In v- dimensional 
space, the Coulomb potential is given by (14.351) . and the corresponding electric field E by 
(I4.34p . with e replaced by e^. The projectile therefore acquires a momentum transfer 



Ap = e p / dtE(h + v p t) , (4.37) 

J — oo 

and it suffers a corresponding change in energy 

AE = . (4.38) 

2m p v ; 

The component of the electric field along the direction of motion v p integrates to zero in 
(I4.37p . while the component normal to the trajectory gives the impulse 

Ap = e p dtE ± {b + v p t)b = ^ / dtb{b 2 + v 2 p t 2 )-»/ 2 b (4.39) 

J —OO ilv— 1 J —OO 

1 b . (4.40) 



ku-2 



14 In this simple example, I am not concerned with hyperbolic orbit corrections and the like; but rather, I 
am tracing the origin of the logarithmic divergence of the Coulomb potential and the pole at v = 3. For 
this, we can work with a hot dilute plasma where a linear trajectory will suffice. 
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The temperature equilibration rate contains a factor involving the cross-section weighted 
energy transfer, and we see that the energy exchange between p and a can be written 

dE f* /' femax db 

— daAE~ — 

dt I L . b v 



doAE~ I (4.41) 



where we have used the fact that da ■ AE ~ b v ~ 2 db ■ /3~ 2 ( y ~ 2 ) ~ db/b u ~ 2 . I will elaborate 
further on this example in the next lecture, but for now one should simply note that the rate 
f)4.4ip implies that large v is dominated by short-distance physics and small v is dominated by 
long-distance physics. Moreover, expression (I4.4ip gives v = 3 as the dividing line between 
these two regions. To see this, note that for v > 3 the impact parameter integral is not 
sensitive to the large distance cut off, and we may simply take the limit 6 max — > oo to obtain 
f°° db b 3 ~ u 

v > 3 : / — - = => UV dominant and pole at v = 3 . (4.42) 

Conversely, for v < 3, we can set b min = 0, with 

v < 3 : / —r = -^5*. =>. IR dominant and pole at v = 3 . (4.43) 



o 



b u ~ 2 3 -is 

The results displayed are the dominant forms in the two regions of spatial dimensionality 
about v = 3. In either case, the stopping power contains a pole 1/(^ — 3). 

There are a number of important consequences arising from the UV and IR behaviors of 
the Coulomb potential (14.351) . I will discuss this more fully in Section I VI Al but for now recall 



that the derivation of the Boltzmann equation, as presented in Ref. 14J for example, breaks 
down for the Coulomb potential in three spatial dimensions because of the aforementioned 
infrared singularity. However, in v > 3 the "textbook derivation" of the Boltzmann equation 
with the Coulomb potential (14.351) is finite and completely rigorous. The simple pole at v — 3 
in the scattering kernel corresponds to an IR divergence because of the long-range nature 
of the Coulomb force in three dimensions. Furthermore, because dimensions greater than 
three enhance the UV physics, the classical Born, Bogoliubov, Green, Kirkwood, and Yvon 
(BBGKY) hierarchy reduces to the Boltzmann equation to leading order in g (or to leading 
order in the number density) when v > 3. A similar reduction from BBGKY holds for the 
Lenard-Balescu equation in v < 3, and the "textbook derivation" |l5| is also rigorous in 
these dimensions. In v = 3, the derivations of the Boltzmann and Lenard-Balescu equations 
break down for the Coulomb potential. This is not because the v = 3 version of the BBGKY 
hierarchy is divergent, but because the truncation procedure that leads to the Boltzmann 
and Lenard-Balescu equations breaks down for the three dimensional Coulomb potential. 15 
Indeed, BBGKY is completely finite for the Coulomb potential in v = 3, albeit completely 



useless for our purposes. Section II of BPS [1] provides more details, especially the two 
paragraphs between Eqs. (2.7) and (2.8). 



15 



For example, in deriving the Boltzmann equation from BBGKY, one invokes a principle of uncorrelated 
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V. THE LAMB SHIFT AND THE COULOMB LOGARITHM 



The Lamb shift is interesting for us because it provides another connection with the 
Coulomb logarithm, both in the historical details and in much of the physics. The Lamb 
shift is a small energy split in the otherwise degenerate 2S and 2P states of total angular 
momentum j = 1/2. The measured value is about APi amb — 4.4 x 10~ 6 eV, with the 2S , 1 / 2 
state lying above the 2P\/2 state. Therefore, when an electron makes a transition from 
the 5-state to the P-state, it emits a microwave photon of frequency Az/i am b — 1060 MHz. 
Calculating the observed value of the Lamb shift was the first great success of quantum 
electrodynamics (QED), the relativistic quantum theory of light and matter. 16 

Dirac's relativistic model of the hydrogen atom, i.e. his relativistic theory of the electron 
in a Coulomb potential, made the prediction that the energy levels of the hydrogen atom 
(neglecting the hyperfine structure) depend only upon the principal quantum number n and 
the total angular momentum j (the sum of the orbital momentum I and the spin s = 1/2 
of the electron). In particular, Dirac predicted that the 2Si/ 2 and 2Pi/ 2 states should be 
degenerate. At the 1947 Shelter Island conference in New York, W. Lamb and R. Retherford 
announced the results of their highly sensitive experiment measuring the emission frequency 
of photons in a 2S-2P transition, thereby establishing the unequivocal experimental existence 
of the Lamb shift. 17 While almost degenerate compared to the binding energy of the hydrogen 
atom, today's accepted experimental splitting is (l6j 

AP lamb = E 2Sl/2 - E 2Pl/2 = 4.374898(7) x 10~ 6 eV (5.1) 
A E 

Az/ lamb ee — ^ = 1057.845(9) MHz , (5.2) 

where Plank's constant in the form h = 4.135 66733(10) x 10~ 15 eV-s is the conversion factor 
between energy and frequency. Hans Bethe had attended the Shelter Island conference, and 



collisions to replace the two-point correlation function by the product of two one-point functions, and this 
leads to the infrared divergence. This means, of course, that it is the long-distance correlations themselves 
that tame the IR divergence of the Boltzmann equation, correlations that are neglected by the truncation 
process. As a related point, this procedure imposes an implicit initial condition, thereby providing a 
direction of time for the Boltzmann equation (and this occurs even for short-range potentials). 

16 The other early success of QED, which followed soon after the Lamb shift, was calculating the magnetic 
dipole moment of the electron. It is customary to write the magnetic moment of the electron \x e in terms of 
the Bohr magneton /i B = eh/2m e by introducing the dimcnsionlcss g-factor: /i e = g e /i B . Using relativistic 
single-particle quantum mechanics, in 1928 Dirac predicted g e = 2 (exactly). In 1948 Schwinger used QED 
to calculate the radiative corrections, and he found g e /2 = 1 + a/2n + 0(a 2 ) = 1.0011614, which was in 
excellent agreement with experiment. Today, the electron's magnetic dipole moment has been calculated 
to include a terms, and is in agreement with experiment to 10 significant figures, the most accurately 
verified quantity in the history of physics. 

17 In the 1930s, S. Pasternak analyzed experimental data suggesting that such an energy split might exist; 
however, the systematic error was as large as the energy splitting itself. 
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on the train ride back he performed the first calculation of the Lamb shift, finding the value 
A^bethe = 1040 MHz. Bethe's calculation was recognized as being only a rough approxima- 
tion, neglecting high energy relativistic effects, but its close agreement with experiment was 
cause for optimism. 

The dominant contribution to the Lamb shift comes from the bound 25* and 2P electrons 
exchanging virtual photons with the atomic nucleus. 18 These radiative corrections effectively 
smear the point-like nature of the nucleus, thereby altering the atomic energy levels (the elec- 
tron no longer "sees" a pure Coulomb potential). Any possible number of photon exchanges 
with the nucleus are permitted, from a single high-energy photon (hard/UV physics) to many 
low-energy photons (soft/IR physics), and this means there are two disparate but competing 
energy scales in the problem. In a manner similar to the Coulomb energy- loss rate (12 .11) . the 
Lamb shift takes the form 

AE lamh = K\J^\ . (5.3) 



The UV scale is set by -E max ~ m e c 2 = 511 keV and the IR scale E min is set by the binding 
energy of the hydrogen atom 

E ° = \ (i^) me = l a2 me ° 2 = 13 - 6 eV • (5 - 4) 

To understand the origin of these scales, note that -E max takes its value from the energy 
at which relativistic effects for the electron become important (the rest-mass energy of the 
electron), while E min takes its value from the only low-energy scale in the problem (namely, 
the binding energy of the atom). The soft photon IR interactions can be handled by nonrela- 
tivistic means, while the hard UV photons require a more complicated relativistic treatment. 
For this reason, Bethe used the simpler nonrelativistic formalism, cutting his calculation off 
at the relativistic energy scale E max ~ m e c 2 at which his formalism broke down. This is 
akin to cutting the calculation of the energy exchange rate off at some small distance scale 
b min . Bethe then concentrated on the low energy theory, which is analogous to looking at 
the Lenard-Balescu equation in the rate calculation. Bethe found a logarithmic UV diver- 
gence, in the same way the Lenard-Balescu equation has a UV divergence, but he was able 
to regularize the infinity by applying a technique known as mass renormalization, thereby 
rendering his calculation finite. His calculation was completely rigorous at low energies, and 



Bethe found the coefficient K and the low energy cutoff -Emm, both exactly |l7j ]: 



a 3 

K=—E Q and E min = A E with £ max ~m e c . (5.5) 

37T 



18 There are also several other subdominant mechanisms, such as the vacuum polarization of the photon 
and the anomalous magnetic dipole moment of the electron. These give, respectively, the contributions 
A< a ^ b = -27 MHz and Az^fb = +68 MHz. 
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The binding energy Eq of the hydrogen atom is given by (15.41) . and while I will not write 
it down, the coefficient A Q is rigorously defined in terms of sums of matrix elements of the 
various intermediate states. These matrix elements are sufficiently complicated that one can 
only calculate them numerically, and to three significant figures Bethe found A$ = 17.8. 
Bethe's exact calculation of E m i n is akin to an exact calculation of 6 max in the rate problem. 
This, however, is a point where the analogy is not precise: Bethe was able to exactly calculate 
E min through the regularization procedure of mass renormalization, which has no counterpart 
in plasma physics. On the other hand, Lyman Spitzer was only able to estimate the value 
of the maximum impact parameter in the Coulomb logarithm to be of order 6 max ~ k,~ x . 
Bethe was able to perform this feat because QED is a fundamental theory of nature, while 
the Boltzmann equation is not. Using (15.41) . we can express (15.31) as 



where a = e 2 /4nhc = 1/137.036 is the fine structure constant. Bethe's calculation therefore 
gives A£f amb = 4.3 x 10" 6 eV or Az/* mb = 1040 MHz. 

Bethe did not, however, calculate the exact coefficient under the logarithm. Since his cal- 
culation broke down at relativistic energies, he used the somewhat ad hoc value -E max = w e c 2 
(rather than some multiple thereof) for the UV cutoff. It turns out that the high energy 
corrections to (15.61) are rather small, so Bethe's result was perhaps more accurate than war- 
ranted. Since Bethe only estimated the maximum energy cutoff, his result ( I5.6P is only 
accurate to leading order in a, and the constant under the logarithm required a more accu- 
rate treatment. The analogy further continues: Bethe calculated the Lamb shift to leading 
order, just as Spitzer calculated the Coulomb logarithm to leading order. Both calculations 
found the correct coefficient in front of the logarithm, and the order of magnitude of the 
argument of the logarithm. These calculations failed to extract the exact constant under the 
logarithm, although Bethe managed to find the exact expression for the low-energy cutoff. 
In this sense, Bethe is the Spitzer of the Lamb shift. 

Shortly after this, R. Feynman and J. Schwinger independently calculated the high energy 
contribution using their respective relativistically covariant formalisms, 19 and their calcula- 
tions were exact to leading and next-to-leading order in a. After adding their contribution 
of the high energy corrections to Bethe's low energy form, the calculation of the Lamb shift 
was complete and the constant under the logarithm was fully determined. Or so it would 

19 F. Dyson soon proved that Schwinger's mathematically rigorous formalism was equivalent to Feynman's 
more intuitive but easier to use approach. Along with the Japanese physicist S. Tomonaga, Schwinger 
and Feynman shared the 1965 Nobel prize "for their fundamental work in quantum electrodynamics, with 
deep-ploughing consequences for the physics of elementary particles." 



lamb 




(5.6) 
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seem. Simultaneously, J. French and V. Weisskopf completed an independent calculation of 
the Lamb shift that disagreed with Feynman-Schwinger, albeit only slightly in the coefficient 
under the logarithm. The plot thickened as yet another calculation of the Lamb shift was 
completed that agreed with French- Weisskopf, this time from N. Kroll and Lamb himself. 20 
So now we seem to have a real problem: the Schwinger-Feynman calculations agree with 
each other, but disagree with the French- Weisskopf and Kroll-Lamb calculations. Schwinger 
and Feynman had developed independent but equivalent formalisms that manifestly exhib- 
ited the relativistic covariance of the theory, while the other four had used a cumbersome 
formalism developed in the 1930's (now called old-fashioned perturbation theory). As it turns 
out, Schwinger and Feynman had made the same subtle mistake, and Weisskopf 's calculation 
was correct. The contribution to the Lamb shift we have been considering has the form 

a 3 /, fm e c 2 l 91\ , . 

where En = 17.8£7o is Bethe's low energy result ( 15. 51) . As with the BPS Coulomb logarithm, 
we can bring the additive constant in (15. 7p inside the logarithm, and using (15. 4p we can 
express ( 15 . T[) as 

a 5 

A£i amb = m e c 2 ln{Ca} with C = 2.98 . (5.8) 

3n 

In terms of the Bethe's coefficient C B of (15. 6p . the complete Lamb shift coefficient is 
C = \[2 e~ 91 / 240 C B (apart from the photon vacuum polarization and the electron dipole 
moment corrections, which were mentioned in footnote [TH but will not concern us further). 

French spent the next year tracing down the origins of the discrepancy between the 
calculations. He could not find an error in the Feynman-Schwinger high energy calculation, 
nor was there an error in Bethe's low energy calculation. Instead, the error was in the way 
the high-energy calculation of Feynman-Schwinger was "joined" onto the low-energy result 
of Bethe. The high-energy calculation is also infinite and in need of regularization, but this 
time it contains an IR divergence. This is analogous to the Boltzmann equation containing 
the correct short-distance physics, but nonetheless diverging at long distances. Therefore, 
Feynman-Schwinger had to introduce an intermediate step in which they regulated their high 
energy theories in the infrared. They chose to do this by assigning a small mass to the photon, 
and then taking this mass to zero at the end of the calculation. While this is a common IR 
regularization scheme in QED, it is incommensurate with Bethe's regularization scheme in 
his low energy calculation, where he used a simple cutoff procedure in a high-energy integral. 

20 Lamb was both an experimentalist and a theorist, and he received the 1955 Nobel Prize for his "discoveries 
concerning the fine structure of the hydrogen spectrum." 
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French- Weisskopf and Kroll-Lamb had gotten the correct result because they calculated both 
the low-energy and the high-energy contributions using the same formalism, and therefore 
with consistent regularization schemes for both short and long distances. For more on 
the consequences of regulating large and small scales in an incommensurate manner, see 
Feynamn's footnote 13 in Ref. |18j of this lecture, which I have quoted in the bibliography in 
its entirety. That one must calculate the large and small scales in exactly the same manner 
is not a minor point, as underscored by the stature of the physicists who failed to realize 
its importance when they first calculated the Lamb shift from the then fledgling theory of 
quantum electrodynamics (QED). 

There are a number of consistent regularization schemes in use today in QED and other 
field theories of nature, with dimensional regularization being one of the most popular and 
easy to use. These regularization and renormalization schemes have allowed us to calculate 
a great many experimentally verified quantities, to extremely high precision, and there is no 
longer any doubt in their correctness. While extending the dimension of space to complex 
values might at first seem unsettling, this procedure works. When dimensional continu- 
ation was first introduced into quantum field theory, there were strong reactions against 
it. However, calculations performed with this method agree with calculations using other 
regularization schemes, and more important, with experiment. In time, particle physicists 
learned to accept the notion that one can perform correct three dimensional calculations by 
working in arbitrary complex dimensions. For our concerns, we note that Refs. jsj and 16] 
have indeed calculated the Lamb shift using the method of dimensional continuation. The 
problem of Feynman-Schwinger is avoided, and dimensional continuation gives the correct 
experimentally observed result with much less effort than more traditional methods. As the 
above discussion illustrates, this is no small achievement. 

In summary, dimensional continuation is powerful because (i) it is a consistent regulariza- 
tion scheme that (ii) lends itself to a perturbative analysis, and (iii) requires relatively simple 
(or at least easily learned) calculational tools. De Witt's calculation [7] was certainly consis- 
tent, in that it matched the long and short-distance physics commensurately [to order 0(g 3 )]. 
This is because he starts with a finite regulated theory that treats long- and short-physics 
together within a single framework, albeit with a regularization scheme that does not lend 
itself to a systematic perturbative analysis. By comparison with the BPS result, we know 
DeWitt was accurate to order g 2 , inclusive. But, as illustrated by Feynman's footnote 13 of 



Ref. 18|, and by the story of the Lamb shift above, any attempt at treating the long- and 
short- distance physics by separate regularization schemes (as much of the plasma literature 
currently does), will likely miss the very constants they are trying to calculate. Recall, it 
took two correct calculations to find the error in Feynman's single calculation. Dimensional 
continuation performs all necessary book keeping, at both small and large scales, and in a 
manner that affords simple calculations and perturbative expansions. 
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VI. CALCULATING THE RATE SYSTEMATICALLY WITH DIMENSIONAL 
CONTINUATION 



A. Dimensional Reduction of BBGKY 

Let us return to the rate equation (jl.7p in the light of the apparatus of dimensional 
continuation that we have constructed. Since we are interested in spatially uniform plas- 
mas, we will only consider particle distributions that are functions of the momentum, or 
equivalently the velocity. Let \ v denote a //-dimensional velocity vector with components ve 
for £ = 1, • • ■ , u, and define a //-dimensional distribution function f v {y v , i) by 

d u v f u (v u ,t) = number of particles in a hyper-volume d v v about w u at time t . (6.1) 

Then the generalization of the three dimensional result (11.71) to //-dimensions would be 

d£ v f ju.. 1 _ . 2 d L 



it J d " v 2 m - "5P V "' ( ' ' < 6 ' 2 > 

where the square of the velocity in (16.21) is v 1 = v u ■ v v = Yle=i v i- As previously men- 
tioned, the standard textbook calculation of the Boltzmann equation goes through without 
an infrared divergent scattering kernel when v > 3. For now, I will write this equation in 
schematic form as 

df 

^ + v,-V/„ = B u [f) :u>3. (6.3) 

Note that V is the //-dimensional gradient and v^-V/^ = ^Y^i=\ v i^fvld^i- In dimensions 
v < 3, the Lenard-Balescu equation is ultraviolet finite, and we have 

| + v„-V/ v = L4/] :u<3. (6.4) 

The explicit form of the scattering kernels B v [f] and L v [f] will not be required until the 
next lecture. As we shall see, the scattering kernels B u and L v are simply the obvious gen- 
eralizations of their three dimensional counter parts: momentum and wave-number vectors 
live in //-dimensions, and the scattering is produced by the //-dimensional Coulomb potential 
(14.35}) . Calculations using (16.31) in v > 3 and (16.41) in v < 3, respectively, are completely 
finite. 

In exactly the same manner, one can generalize the BBGKY hierarchy to an arbitrary 
number of dimensions, and this will be the starting point for my treatment of dimensional 
continuation. As we discussed in the closing paragraph of the last section, when the number 
of spatial dimensions is greater than three, BBGKY reduces to the Boltzmann equation (16.31) 
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FIG. 4: For ^>3 the "textbook derivation" of the Boltzmann equation for a Coulomb potential is rigorous; 
furthermore, the BBGKY hierarchy reduces to the Boltzmann equation to leading order in g. A similar 
reduction from the BBGKY hierarchy holds for the Lenard-Balescu equation in v < 3, and the "textbook 
derivation" is also rigorous in these dimensions. In v = 3, the derivations of the Boltzmann and Lenard- 
Balescu equations break down for the Coulomb potential. This is not because the three dimensional BBGKY 
hierarchy breaks down, but because the Boltzmann and Lenard-Balescu equations break down. Indeed, 
BBGKY is completely finite for the Coulomb potential in v = 3, albeit completely useless for our purposes. 

to leading order in the plasma coupling g. Therefore, when v > 3, to leading order in g the 
rate becomes 

^L = Jd»v^mvlB v [f] : u>3. (6.5) 

Conversely, in dimensions v < 3, BBGKY reduces to the Lenard-Balescu equation (16.41) to 
leading order in g, and 

d£< r ... 1 



(Jf d v v -mvi L u [f] : u < 3 . (6.6) 

This is illustrated in Fig. HI The "greater-than" and "less-than" superscripts are to remind 
us that the integrals in the rate calculations (16. 5ft and (16. 6ft are to be preformed in v > 3 and 
v < 3, respectively. Because of spatial uniformity, I have set the convective terms in (16.31) 
and (16.4p to zero, v u - V/^ = 0. Using the appropriate Coulomb potential (I4.35P for V v (r) in 
the scattering kernels of (16.51) and (16.61) . the integrals now converge, and they are calculated 
exactly in Sections 7 and 8 of BPS In the next lecture, we will calculate these integrals 
in the Born approximation. 

Note that this is a first-principles derivation of the rates (16.51) and (16.61) in their respective 
dimensions v. Let me reiterate the argument once again, although with a slightly different 
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emphasis. For simplicity, we consider the purely classical regime first, adding quantum me- 
chanics in a moment. The classical BBGKY hierarchy for the Coulomb potential in three 
spatial dimensions is well defined and finite. We run into trouble only when we attempt 
to truncate the hierarchy and derive lower-order kinetic equations, such as the Boltzmann 
and the Lenard-Balescu equations. The necessity of truncated equations is of course clear 
(even the three-body problem cannot be solved analytically). Unfortunately, however, the 
Coulomb potential in three dimensions produces divergent scattering kernels in these trun- 
cated equations. Rather than creating a model of the ostensibly divergent scattering kernel, 
we shall instead systematically regulate the divergences by letting the spatial dimension 
depart from its empirically measured value of three. Logarithmic divergences in three di- 
mensions then become simple poles of the form — 3) in arbitrary dimensions. As with 
the regularization procedure of quantum field theory, our starting point here is a well-defined 
and finite theory, albeit in v dimensions, regularized in a consistent fashion at all length and 
energy scales. 

Let us now return to classical BBGKY, but this time in a spatial dimension v of arbitrary 
positive integer value (we are not yet considering continuous values of v). At first sight, the 
hierarchy equations in v dimensions are just as useless as those in three dimensions — there 
are simply too many of them to solve. However, if we are willing to work to leading order ac- 
curacy in the plasma coupling g, which is quite accurate for a weakly to moderately coupled 
plasma, then: (i) in v > 3 we can truncate BBGKY to the //-dimensional Boltzmann equa- 
tion, and (ii) in v < 3 we can truncate BBGKY to the //-dimensional Lenard-Balescu equa- 
tion. Quantum scattering effects in the plasma will not modify the //-dimensional Lenard- 
Balescu equation, but they will modify the //-dimensional Boltzmann equation. Since we 
require the leading order term in the rate to be exact, we must include two-body quantum 
scattering effects exactly (but no more than two-body effects, since these are subleading 
in g — in fact, three-body and higher correlations and scattering effects, both classical and 
quantum, can and should be neglected to leading order in g). Two-body quantum effects 
can be accounted for by replacing the classical cross section in the Boltzmann scattering 
kernel with the the two-body quantum cross section (this can be performed exactly since all 
scattering phase shifts 5e are known for the Coulomb potential in three dimensions). These 
ideas will be illustrated in complete detail in the next lecture for a particularly useful but 
simple case. 
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B. Calculating the Leading Order Term 

We now return specifically to the electron-ion temperature equilibration rate of (jl.2p . To 
obtain the leading order in g behavior when v > 3, we calculate the rate using (16.51) . As 
our calculation in the next lecture will reveal, this rate is proportional to g 2 (or the number 
density) and takes the form 21 

-f - = H(u; r]) + Q(u - 3) : LO in g when v > 3 , (6.7) 

where, for defmiteness, we have restored the electron- ion subscript to the rate, as in (11. ip . 
In the next lecture we shall calculate H(u; rj) explicitly, but for now it will suffice to note 
that H depends upon the spatial dimension v and the quantum parameter 77. The simple 
pole at v — 3 reflects the long-distance or infrared divergence of the Coulomb potential in 
three spatial dimensions, and it arises from an integral over the radial coordinate of the 
//-dimensional Coulomb potential (I4.35p . In a similar manner, the leading order behavior in 
dimensions v < 3 is given by (16.61) . In the next lecture, we will see that this takes the form 

f]£ < a v ~ x 

_« = G{v) f + 0(3 -v) : LO in g when v < 3 . (6.8) 

There is no //-dependence in G{y) since the leading order long-distance physics is purely 
classical. In the next lecture, we will preform the integrals to establish (16. 7p and ( 16.8ft . 
thereby calculating the coefficients H{y) and G{y) exactly (for notational simplicity, I will 
hereafter drop the ^-dependence from H). Note that (16. 8p also contains a simple pole at 
v = 3 arising from an integration of ( 14.35ft . but this time the pole corresponds to missing 
short- distance physics or the ultraviolet divergence of the Lenard-Balescu equation in three 
dimensions. In general, if a three-dimensional integral diverges, then the corresponding 
v- dimensional integral will typically contain a simple pole of the form \j{v — 3). In this way 
we can transform a divergent integral into a convergent quantity that we can manipulate, 
and this is how we shall regularize the divergent Boltzmann and Lenard-Balescu equations 
in three dimensions. 



As we have already pointed out in footnote [U the parameter g = e 2 K D /47rT given by (|4.2|) is dimensionless 
only for v — 3. From definition (|4.ip and the Coulomb potential (|4.35j) . the proper dimensionless expansion 
parameter should be g v — e 2 K^~ 2 /C V T. However, since g u — g ■ '(47T / 'C '„) , we can absorb factors 
of Rq~ 3 (47t/C„) into any accompanying coefficient, such as H(i/;r]) in (|6.7[) . Therefore, powers of the 
dimensionless coupling g v also count powers of the three dimensional coupling g. We are therefore free to 
think of the expansion, even in v dimensions, in terms of the three dimensional parameter g. 
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analytically continue 
around the v = 3 pole 




d£</dt _ d£</dt 



g 2-(3-u) 3 g 2-(3-u) = g 2+(u-3) 



4 ^ 

„2— 3| „2+|i;-3| 

LO: large when g <C 1 NLO: small when <? ^ 1 

FIG. 5: The analytic continuation of dS^/dt from v < 3 to the region ^ > 3: the same expression can be 
used for d£</dt throughout the complex plane since the pole at v — 3 can easily be avoided. Note that the 
quantity dE^/dt ~ g 2 +( u ~ 3 ) is leading order in g for v < 3. However, upon analytically continuing to v > 3 
we find that dS^/dt ~ g 2+ ' ly_3 ', which is ncxt-to-lcading order in 5 relative to d£>/dt ~ g 2 . 



C. Next-to-Leading Order from Leading Order via Analytic Continuation 

Since the rates d£>/dt of (16.71) and d£</dt of (I6.8j) were calculated in mutually exclusive 
dimensional regimes, one might think that they cannot be compared. However, even though 
( 16. 8p was originally calculated in v < 3 for integer values of u, we can analytically continue 
the result to complex values of v (in the same way that the factorial function, which operates 
on positive integers, can be generalized to the Gamma function, which operates on the whole 
complex plane). In fact, if we continue d£</dt to real values of v with v > 3, then we can 
directly compare (16 .7h and ( I6.8p . Upon writing the (/-dependence of ( 16. 8 h as g 2+ ( u ~ 3 \ when 
v > 3 we see that ( 16. 8j) is higher order in g than ( 16. 7\} . By power counting arguments, there 
are no powers of g in ( 16.71) between g 2 and g v ~ l for v > 3, and therefore (16.81) indeed provides 
the correct next-to-leading order term in g when the dimension is analytically continued to 
v > 3, 

d£< 2+(^-3) 

_« = -G{v) — + 0(u - 3) : NLO in g when v > 3 , (6.9) 

where I have written the exponent of g as 2 + {y— 3) rather than v — \. This is illustrated 
in Fig. [51 

As we shall see in the next section, the individual pole-terms in ( 16. 7p and (16.91) will cancel 
giving a finite result when the leading and next-to-leading order terms are added. The 
resulting finite quantity will therefore be accurate to leading and next-to-leading order in 
g as the v —>■ 3 limit is taken. Alternately, we could have analytically continued (16. 7p to 
dimensions v < 3, where it would become next-to- leading order relative to (16.81) . In either 
case, the leading and next-to-leading order contribution can be found by simply adding (16.71) 
and (16.81) and then taking v — > 3. 
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D. Returning to Three Dimensions 



To find the three dimensional rate d£ el /dt, accurate to leading and next-to-leading order 
in the plasma coupling g, we add the leading order expression (16. 7p for d£> 1 /dt to the next- 
to- leading order expression (16.91) for d£</dt, and then take the limit v — > 3 + : 



d£ a ,. 
= hm 

at i>->3+ 



d£ e i _|_ ^el 

dt dt 



+ 0(g 3 ) . (6.10) 



For the same reasons as given in Section IIIIEl this does not lead to any form of "double 
counting." Instead, we are simply adding the next-to-leading order term (I6.9P to the leading 
order term (16.71) at a common value of v > 3. 

Since (I6.10p lies at the heart of dimensional continuation, allow me review the reasoning 
behind this expression one final time. Recall from (16. 7p and (I6.8P that d£>/dt and d£</dt 
are both leading order in g for v > 3 and v < 3, respectively. Since these functions were 
calculated for mutually exclusive values of u, they must be analytically continued to the 
same value of v for purposes of comparison. In equation (I6.10p . we have chosen to continue 
d£<Jdt to v > 3, which takes the same algebraic form as it did for v < 3. The analytic 
continuation has the effect of rendering d£</dt subleading in g relative to d£>/dt. There 
are no powers of g between these two terms for any value of u, and therefore the limiting 
procedure v — ► 3 + gives the three dimensional rate exactly to leading and next-to-leading 
order accuracy in g. We have therefore found the leading order and next-to- leading order in 
g contributions for v > 3, 



If + H = ^3 [ H{V) ^ ~ G{U) g2+iV ~ 3) ] + ° {V ~ 3) ' (6 " n) 

I have not indicated the higher order error in g, but instead I am emphasizing here the 
error associated with the //-expansion. That is to say, ( 16 .111) contains an implicit error that 
approaches 0(g 3 ) as v — > 3 + , while I have chosen to explicitly display the first-order error 
in v — 3. In the sum (16. lip , the error in g remains nonzero in the three dimensional limit, 
while the z/-error in (16. lip vanishes as v — > 3. 

In taking the v — > 3 limit, we must keep in mind that terms proportional to v — 3 will 
give a non- vanishing result when multiplied by the pole l/{y — 3). Upon expanding the 
g-dependence to linear order in e = v— 3 we find 

g e = exp{ln(# e )} = exp{eln#} = 1 + e In # + 0(e) 2 , (6.12) 

which can be written as 

j = i + ln# + 0(e) . (6.13) 
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We must now expand the coefficients H{y) and G{v) in powers of e = v — 3. As we shall 
see in the next lecture, and this is a crucial point, the leading order in e terms are identical, 
and so the expansions take the form 

H{v) = -A + e^ + Oie 2 ) (6.14) 
G(u) = -A + ed + 0(e 2 ) . (6.15) 

In the next lecture, we will perform the integrals in (16.51) and (16. 6p . thereby allowing us 
to exactly compute A = H(u = 3) = G{y = 3) and the linear terms Hi = H'(u = 3) and 
G\ = G'(u = 3). The remaining procedure is now straightforward. Upon using (I6.13I) - (I6.15I) 
in (16.71) and (16. 9p . the leading and next-to-leading order terms in g now take the form 

d -§- = 9 2 + H l9 ' + 0(u - 3; g*) (6.16) 

= ^g 2 -AgHng-G ig 2 + 0(v-3;gZ). (6.17) 

I have now indicated here the error associated with the ^-expansion. Because the leading 
order terms in (16.141) and (16.151) are equal, the simple poles in (16.161) and (16.171) cancel. 
Therefore, the limit v — > 3 + of expression (16.111) gives the rate 

^ = -AgHng + Bg> + 0(g*), (6.18) 

with B = H i — G 2 , in agreement with (14.71) . In this way, BPS has calculated the energy 
exchange accurate to leading order and next-to- leading order in g. 



VII. SOME CLOSING REMARKS 



A. Summary 

I hope this account of the simple field theoretic apparatus necessary for a reading of BPS 
has been helpful. Having understood the reasoning behind BPS, especially the manner by 
which analytic continuation turns a leading order result into a next-to-leading order result, 
we can summarize the BPS procedure by the following prescription: 

1. Calculate the rate d£ > /dt in the regime v > 3 using the Boltzmann equation, including 
quantum corrections to all orders in rj. This captures the leading order physics in 
dimensions greater than three. 

2. Calculate the corresponding rate dS < /dt in the regime v < 3 using the Lenard-Balescu 
equation. This physics is classical, and captures the leading order behavior in dimen- 
sions less than three. 
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3. Add the terms dS > /dt and d£ < /dt and take the v — > 3 limit. This gives the energy 
exchange rate in three dimensions accurate to leading and next-to- leading order in g, 
dS/dt = — Ag 2 \n{C(r])g} + 0(g 3 ). That is to say, this gives the coefficients A and C 
exactly. 

This prescription provides complete analytic expressions for the coefficients A and C(rj), 
with C(r]) being of particular interest. In the language of the Coulomb logarithm, we write 
mA coul O,77) = - \n{C(rj)g}. 



B. Further Context 

The fundamental interactions of nature can be expressed as quantum field theories (except 
for gravity 22 , which to date has only been expressed as a classical field theory, i.e. general 
relativity). For this reason, quantum field theory reveals something very deep about the 
structure of nature. However, for our purposes, we only need to think of quantum field 
theory as an elegant collection of tools for packaging and solving many-body problems, 
and plasma physics concerns itself with the many-body problem par excellence. It should 
therefore come as little surprise that quantum field theory can be useful for plasma physics. 
Reference 
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is a nice example of cross fertilization between two quite different branches 
of physics, and it is gratifying for someone who has worked in both particle and plasma 
physics that such seemingly different subjects can inform one another. Perhaps the most 
dramatic difference between particle and plasma physics lies with their respective methods 
and outlooks, and not so much the subject matter itself. A number of plasma physicist 
have explained to me that that plasma physics is not an exact science. While the nature of 
the subject might render a first-principles approach limited, particularly for subtopics like 
magneto-hydrodynamics in tokamaks or self-organized behavior in strongly coupled plasmas, 
I believe the plasma physicist could still benefit from the more rigorous outlook of particle 
physics. Conversely, the particle physicist could benefit from plasma physics. In fact, until 
new high energy experiments come on-line, plasma physics may have more to offer particle 
physics than the reverse. 23 Until recently, high energy physics was a spectacular field in 



22 There have of course been numerous attempts at quantizing gravity, but since none of these has yet 
produced a full fledged theory of quantum gravity, I do not count them here. In fact, it might not be 
possible to describe gravity with a field theory on the quantum level. This would be true, for example, if 
string theory were to provide the theory of quantum gravity, since a string theory is qualitatively different 
from a field theory. 

23 Particle physics is currently in crisis brought on from a dearth, and indeed a complete absence, of new 
experimental results. The only hope for high energy physics seems to lie with the Large Hadron Collider 
at CERN, or its possible successors. Compare this situation with the current Renaissance in observational 
cosmology and astrophysics, which has more and better quality data than at any time during its history. 
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which to work, but, with a few exceptions such as lattice gauge theory and RHIC physics 
at BNL, the field seems to have been milked dry. High energy physics is in dire need of 
direct experimental input. Plasma physics, on the other hand, still has many interesting 
and quite challenging problems of experimental and astrophysical relevance, and a number 
of these seem tailor-made for the field theorist. The particle physicist could therefore benefit 
by an excursion (or even a longer stay) in plasma physics. Hence, my hidden agenda in 
these lectures: it would be quite nice to entice a few particle theorists to work on some of 
the interesting problems of plasma physics. It would also be nice to convince a few plasma 
physicists that particle theory has something to offer them as well. 
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That the result given in B [R.P. Feynman, Phys. Rev. 74, 1430 (1948)] in Eq. (19) was in error 
was repeatedly pointed out to the author, in private communication, by V.F. Weisskopf and 
J.B. French, as their calculation completed simultaneously with the author's early in 1948 gave 
a different result. French has finally shown that although the expression for the radiationless 
scattering in B, Eq. (18) or (24) above is correct, it was incorrectly joined onto Bethe's non- 
relativistic result. He shows that the relation ln2A; max — 1 = lnA m i n used by the author should 
have been ln2fe max — 5/6 = lnA m ; n . This results in adding a term —(1/6) to the logarithm in B, 
Eq. (19) so that the result now agrees with that of J.B. French and VW. Weisskopf, Phys. 
Rev. 75, 1240 (1949) and N.H. Kroll and W.E. Lamb, Phys. Rev. 75, 388 (1949). The author 
feels unhappily responsible for the very considerable delay in the publication of French's result 
occasioned by this error. This footnote is appropriately numbered. 
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